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ABSTRACT 

We calculate the X-ray emission from both constant and time evolving 
shocked fast winds blown by the central stars of planetary nebulae (PNs) and 
compare with observations. Using spherically symmetric numerical simulations 
with radiative cooling, we calculate the flow structure, and the X-ray tempera- 
ture and luminosity of the hot bubble formed by the shocked fast wind. We find 
that a constant fast wind gives results that are very close to those obtained from 
the self-similar solution. We show that in order for a fast shocked wind to explain 
the observed X-ray properties of PNs, rapid evolution of the wind is essential. 
More specifically, the mass loss rate of the fast wind should be high early on 
when the speed is ~ 300 — 700 km s _1 , and then it needs to drop drastically by 
the time the PN age reaches ~ 1000 yr. This implies that the central star has a 
very short pre-PN (post-AGB) phase. 

Subject headings: Subject headings: stars: mass loss stars: winds, outflows planetary 
nebulae: X-ray X-rays: ISM 



1. INTRODUCTION 

In the last seven years people have been trying to understand the nature and the origin 
of the extended, spatially resolved X-ray emission detected in planetary nebulae (PNs) by 
the ChandraX-iay Observatory (CXO), e.g., BD +30°3639 (Kastner et al. 2000; Arnaud et 
al. 1996 detected X-rays in this PN with ASCA), NGC 7027, (Kastner, et al. 2001), NGC 
6543 (Chu et al. 2001), Henize 3-1475 (Sahai et al. 2003), Menzel 3 (Kastner et al. 2003), 
and by the XMM-Newton X-ray Telescope, e.g., NGC 7009 (Guerrero et al. 2002), NGC 
2392 (Guerrero et al. 2005), and MGC 7026 (Gruendl et al. 2004). The X-ray emitting gas 
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can come from shocked fast wind segments that were expelled by the central star during the 
late post-asymptotic giant branch (AGB) phase and/or early PN phase, and/or the X-ray 
emitting gas may result from a collimated fast wind (CFW), or jets, blown in conjunction 
with the companion to the central star during the late AGB phase or early post- AGB phase 
(Soker & Kastner 2003). These interactions play a significant role in shaping PNs (e.g., 
Balick & Frank 2002, and references therein). Therefore, understanding the X-ray emission 
can teach us about the properties of the central fast wind and CFW, and by that shed light 
on the shaping mechanism of PNs. 

One of the major properties of the X-ray emitting gas to be understood is its relatively 
low temperature of ~ 1 — 3 x 10 6 K. The observed velocity of fast winds in PNs is Vf > 
1000 km s^ 1 , implying a post shock temperature of T > 10 7 K, much higher than observed 
in the X-ray emitting gas. Two solutions are possible: (1) The hot (T > 10 7 K) post-shock 
gas is cooled via heat conduction by interaction with the cold, slow-wind material, which is 
at a temperature of T ~ 10 4 K and most X-ray emission comes from the conduction front 
(Soker 1994; Zhekov & Perinotto 1996; Steffen et al. 2005). A similar effect can result 
from mixing of hot post-shock fast- wind gas and cool slow-wind gas (Chu et al. 1997). 
Instability modes near the contact discontinuity can lead to significant mixing of cold and 
hot gas (Stute & Sahai 2006, hereafter SS06). (2) The X-ray emitting gas comes mainly 
from a slower moderate- velocity wind of vj ~ 500 km s" 1 . This flow can be a post- AGB 
wind from the central star (Soker & Kastner 2003; Akashi et al. 2006 hereafter ASB06; 
and SS06), or two opposite jets (or CFW, Soker & Kastner 2003). The idea of a post- AGB 
wind is supported by the analytical calculations of a spherically- symmetric fast wind done 
by ASB06 and based on the self-similar solution of Chevalier & Imamura (1983). The self 
similar solutions, however, are limited in their ability to account for temporal evolution and 
to properly treat radiative cooling. 

In the present paper, we carry out a numerical study for the expected contribution from 
the fast spherical wind blown by the central star. We run numerical simulations solving the 
hydrodynamical equations with a full treatment of radiative cooling, and by setting proper 
initial and boundary conditions of the interacting winds, e.g., a fast wind evolving with time. 
Some of the results obtained here are similar to results obtained very recently by SS06. The 
remainder of the paper is outlined as follows: The numerical method is described in §2. In 
§3 we compare the spherical symmetrical numerical simulations with the self-similar results 
of ASB06. In §4 we present several cases of a time evolving fast wind. A short summary is 
given in §5. 
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2. NUMERICAL METHOD 

The simulations were performed using Virginia Hydrodynamics- 1 (VH-1), a high reso- 
lution multidimensional astrophysical hydrodynamics code developed by John Blondin and 
co-workers (Blondin et al. 1990 ; Stevens et al. 1992; Blondin 1994). The code uses finite- 
difference techniques to solve the equations of an ideal inviscid compressible fluid flow. We 
use 1024 grid points to resolve the calculated region. The distance between adjacent grid 
points increase as the flow expands. Using more or less grid points (e.g., 512) changes the 
results only by a few per cents or less. Radiative cooling was incorporated for gas tempera- 
tures above 10 4 K, using the cooling function A(T) for solar abundances from Sutherland & 
Dopita (1993; their table 6). Radiative cooling should be treated carefully near the contact 
discontinuity. The hot bubble cools slowly because of its low density, while the dense shell 
cools slowly because of its low temperature (T = 10 4 K) where the cooling rate is low. The 
one or two grid points at the interface between the hot bubble and cold shell could share 
high temperature from the bubble and high density from the cold shell. Therefore, radiative 
cooling at the interface may be overestimated. Actually, low resolution in the code might 
in a sense mimic heat conduction. Our code includes a check of this numerical problem. It 
takes for the cooling rate at the contact discontinuity the minimum cooling rate of the two 
zones (either from the cold shell or from the bubble). 

To mimic the ionizing radiation of the central star, we did not let the gas to cool to 
temperatures below 10 4 K when radiative cooling was included. For the initial conditions 
we take a spherically symmetric slow wind with a constant mass loss rate M\ and a constant 
velocity V\ = 10 km s" 1 . The slow wind fills the space around the center from a minimum 
initial radius R c q = R c (t = 0). This implies that the initial density is p\ = M 1 /47Tf 1 r 2 for 
r > R c0 and zero inside this radius. At time t = a spherically symmetric fast wind with a 
mass loss rate of M 2 {t) and a velocity v^t) is turned on close to the center. 

The reasons for setting up a vacuum in the center of the slow wind before the fast wind 
is turned on are: (1) There is a large uncertainty as to the evolution of the wind in the 
transition from the AGB phase to the PN phase. Therefore, there is no clear parameters to 
use for the wind at this stage. (2) This wind period has relatively low mass, as mass loss rate 
is lower than in the AGB phase. On the other hand, it is dense and slow enough to cool very 
rapidly. So practically, it will form a shell inward to the slow wind. As far as the numerical 
procedure is concerned, we can mimic this intermediate wind segment by changing the value 
of Rcq. As we show here, this does not change our results much. 

The results of one run with a non-evolving fast wind of M2(t) = 1.4 x 1O _7 M yr _1 and 
V2{t) = 500 km s" 1 , but without radiative cooling, is shown in Fig. [T] at time t = 400 yr. 
Rco = 3 x 10 15 cm in this run. In all runs t = is when the fast wind starts. In 
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Fig. 1. — Total number density (upper panel) and temperature (lower panel) at t — 400 yr 
from a spherical symmetric constant-flow simulation with no radiative cooling included, and 
with the following parameters: Mass loss rate and velocity of the slow wind: M\{t) = 
7 x 10~ 6 M Q yr^ 1 and Vi(t) = 10 km s -1 . Mass loss rate and velocity of the fast wind: 
M 2 (t) = 1.4xlO -7 M yr -1 andt> 2 (t) = 500 km s -1 . The marked regions are: (a) undisturbed 
central star fast wind; (b) Hot bubble; (c) Shocked slow wind; and (d) Unshocked nebular 
gas (AGB wind). v c and s 2 represent the velocities of the contact discontinuity and inner 
shock, respectively. 
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this run, M\ = 7 x 10~ 6 M Q yr -1 (more examples can be found in SS06.) The fast wind 
collides with the slow wind and two spherical shock waves are formed. One shock moves 
outwards at radius Ri(t) into the undisturbed slow wind, while an inner shock runs into 
the expanding undisturbed fast wind at radius i^W- Between the two shock fronts lies the 
contact discontinuity at radius R c (t) (see also Fig. 2 by Volk & Kwok 1985). The values of 
the density and the temperature jump across the contact discontinuity, while the velocity 
and pressure vary continuously across it. The volume inside the dense shell R2 < r < R c 
(region b) in runs without radiative cooling is filled with hot shocked fast wind. When 
radiative cooling is included, the relatively dense region behind the contact discontinuity 
cools radiatively and forms a cold layer at T ~ 10 4 K of previously shocked fast wind 
material. The volume filled with hot X-ray emitting gas is termed the hot bubble (region 
b). The shocked slow wind in the shell R c < r < R 1 (region c) cools much faster due to its 
high density maintaining a temperature of ~ 10 4 K. 

We are interested in the X-ray luminosity L x and temperature of the hot gas. The 
luminosity is calculated in the energy band 0.2 — 10 keV, as described in ASB06. This 
range is chosen to reflect approximately the sensitivity regime of the Chandra and XMM- 
Newton telescopes and science instruments. As discussed by SS06, for gas temperatures of 
T x < 3 x 10 6 K most of the cooling takes place via emission at lower energies (< 0.2 keV), 
which is outside the range of the instruments we compare our results with. The temperature 
of the X-ray emitting gas in the hot bubble varies with radius (see Fig. [1]). In ASB06 (eq. 
9 therein), we introduced a mean emission- measure weighted temperature T x that would be 
deduced from a single-temperature model for the X-ray spectrum. Here too, we will refer to 
the mean temperature of the hot bubble as T x . 

In the numerical model, the flow relaxes to a more or less steady state after ~ 100 yr. 
For several hundred years, the velocity of the contact discontinuity, v c , and the velocity of 
the inner shock, s 2 > oscillate around their mean values. However, the influence of these 
oscillations on the observed X-ray luminosity and temperature is very small. At the first 
encounter of the fast wind with the slow wind, the temperature of the X-ray emitting gas is 
that of the post-shock gas. After ~ 100 yr, as the post-shock gas expands with the nebula, 
the hot bubble cools, and T x decreases as most of the radiation comes from the region close 
to the contact discontinuity. Some other flow properties at the very early stages of the 
winds interaction process, are discussed by SS06. SS06 compare also the velocity of the 
contact discontinuity in their numerical simulations with analytical expressions. For the run 
presented in Fig. [TJ we find v c = 36.2 km s _1 . This run is referred to as model B5 in ASB06 
and in SS06, with the difference that ASB06 and SS06 took the undisturbed slow wind 
temperature to be very low, T <C 10 4 K, while here T = 10 4 K. The self-similar solution 
with a cold slow wind (T <C 10 4 K) gives v c = 32.6 km s _1 . The higher temperature here 
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(T = 10 4 K) implies more energy in the flow, hence larger v c . When SS06 take T <C 10 4 K 
and include radiative cooling, they find v c = 27 km s -1 , as a consequence of the loss of 
energy by radiation. The actual temperature of the undisturbed slow wind will rise from a 
few 100 K during the pre-PN (PPN) phase to ~ 10 4 K during the PN phase, after ionization 
starts. Neither we nor SS06 treat this full evolution as we do not include photo-heating. 
We prefer to take the temperature as appropriate for the PN phase, while SS06 take the 
PPN-phase temperature. This is not a major issue, as the temperature of the pre-shock 
slow wind affects only slightly the X-ray properties. On another matter, residual differences 
between our results and those of SS06 can be attributed to slight numerical differences in 
fitting the cooling function given by Sutherland & Dopita (1993). 

At early times, the X-ray emission of the PN depends strongly on the slow wind initial 
radius R c q. This is demonstrated in Fig. [21 which shows the PN evolution with time in the 
L x — T x plane for three different initial inner radii R c q. As explained above, the temperature 
is high at very early times, and decreases toward its self similar solution at later times. At 
early times the X-ray luminosity is very low as there is not much gas in the hot bubble. With 
time, the X-ray luminosity increases as more mass is accumulated in the hot bubble (SS06) 
until the effect of expansion takes over, the temperature stabilizes and the X-ray luminosity 
decreases as in the self similar solution (see next section). To reach the self similar behavior 
takes more time for larger values of R c o, but regardless of R c q all simulations converge to 
the same results as seen in Fig. [2j Therefore, we will not discuss the very early times of 
evolution. In most runs, those lasting 5000 yr, we took R c0 = 8 x 10 15 cm. In the shorter 
runs, lasting 2000 yr, we took R c0 = 3 x 10 15 cm. 

In order to demonstrate the effect of radiative cooling, we present in Fig. [3] a comparison 
between simulations with and without radiative cooling. The effect of radiative cooling is 
most pronounced at the densest and coolest shocked fast wind segments, which reside right 
behind the contact discontinuity. As explained above, no gas has a temperature below 10 4 K. 
In the hot bubble, generally 10 5 < T < 10 7 K and the radiative cooling time r coo i is highly 
sensitive to the temperature (~ T 5 / 2 ). This explains the sharp boundary of the radiatively 
cooled region marked as RC in Fig. [3j We note that the region marked RC is found to 
be the point where the radiative cooling time equals six times the age of the flow. Namely, 
regions with cooling time much longer that the flow time had time to cool. The reason why 
these regions had time to cool is that the gas now residing near the contact discontinuity 
(at R c ) was shocked at much smaller radius, and hence its density was higher and radiative 
cooling time much shorter in the past. Therefore, we would expect that the present radiative 
cooling time at R c is longer than the flow age. 
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32.5: 




Fig. 2. — Simulated evolution of PN X-ray emission from t = 10 yr to t = 1000 yr for three 
cases differing by their initial inner radii of the slow wind: R c0 = 3 x 10 15 cm (asterisks); 
Rco = 5 x 10 15 cm (plus signs); R c0 = 1.2 x 10 16 cm (circles). Data points for all runs are 
given at five times as marked on the lower plot. In all runs v\(t) = 10 km s _1 (as in the rest of 
the paper), Mi(t) = 7 x 1O" 6 M yr" 1 , M 2 (t) = 1.4 x 1O~ 7 M yr" 1 , and v 2 (t) = 500 km s -1 . 
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Fig. 3. — Simulation results for density and temperature profiles with and without cooling. 
The solid lines are the same as in Fig. [H except that T = 10 4 K was set as a strict lower 
limit. The dashed lines show the results for the same parameters, but with radiative cooling 
included. The region marked RC is the region where radiative cooling is most effective. It 
is found to be the point where the radiative cooling time equals six times the age of the 
flow. It is seen that radiative cooling reduces the pressure in the hot bubble and hinders its 
expansion. 
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3. CONSTANT FAST WIND 

We run several simulations with a constant, unevolving fast wind, namely, M^it) and 
V2(t) do not vary with time. The results enable a comparison with the self-similar results 
of ASB06 and serve as reference for the evolving fast wind simulations to be discussed in 
the next section. Radiative cooling is included in all runs in this section. We recall that in 
order to account for radiative cooling in the self similar solutions, we simply disregard wind 
segments with cooling times shorter than the PN age (t C ooi < t, see ASB06) in calculating 
the X-ray emission. 

A comparison of the simulation results with the self-similar solutions is presented in 
Fig. H] in terms of the PN evolution in the L x — T x plane. Three cases are shown. The fast 
wind mass loss rate M 2 in units of 1O _7 M yr -1 is marked next to each track inside square 
brackets [ ], while the fast wind speed in units of km s _1 is in parentheses ( ). The mass loss 
rate of the slow wind is Mi = 7 x 10~ 6 M Q yr -1 and the slow wind speed is v± = 10 km s _1 
in all runs. Short straight lines connect the self similar results to the numerical results at 
the same age t. Evidently, the self similar solutions with our approximation for radiative 
cooling provide very good estimates to the X-ray properties of the nebula. In all cases 
presented in Fig. HJ the self similar solutions reproduce T x and L x to within ~ 30% and 
~ 45%, respectively, of the simulation values. For ages of 200 yr and more, the self similar 
results for T x and L x are good to within ~ 12% and ~ 30%. Recall that the simulation 
results at early times are uncertain as discussed in §2. It can be seen in the figure, that 
using less strict criteria for omitting fast-cooling wind segments (e.g., t coo i < 3t) improves 
the late-time self similar results at the price of larger discrepancies at earlier times. For the 
higher velocity simulations (two right cases in Fig. [4]) the radiative cooling removes cooler 
gas more efficiently at early times, hence the temperature is higher than at later times. For 
the slower fast wind case, v 2 = 400 km s _1 , the post shock temperature is low and density 
high, such that radiative cooling is very efficient. Because of that the entire shocked fast 
wind is cooling at early times. This explains its low temperature at early times. 

In all cases, the general behavior of L x oc t~ v with rj ~ 0.9 is obtained. This dependence 
on time can be understood as follows. The velocity of the contact discontinuity v c , is constant, 
hence the volume V, of the hot bubble increases as V = 4:7r(Rl — Rl)/3 ~ 47r(v c t) 3 /3, and 
the density decreases as n oc (M 2 t)/V oc t~ 2 . As the temperature does not change much, the 
emissivity (power per unit volume) e goes as e oc n 2 A(T) oc n 2 oc t -4 , and L x = Ve oc t~ x . 
This simple explanation, and the close fit of the numerical results to the self-similar solution 
with cooling gas removed (Fig. HJ), enhance our confidence in our results. We find R c oc t ' 98 
and the density inside the bubble to vary as p oc t~ 2 , which gives for the mass in the bubble 
M oc t 0,94 . As more mass cools with time, we expect less mass to be in the bubble. SS06 
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find in their numerical simulations r] = 0.37. We have no explanation for this discrepancy, 
but we note that SS06 find the mass in their X-ray bubble to increase as t 15 , rather than 
the simple expectation of ~ t. They explain it by more efficient cooling at early times. 
However, we note that this behavior continues to late times, t > 2000 yr where the contact 
discontinuity is already at r ~ 10 17 cm. 

The results of the numerical simulations with constant fast wind properties strengthen 
our conclusions from self similar calculations (ASB06; see also SS06). Basically, X-ray emis- 
sion from PNs can be accounted for by shocked wind segments that were expelled during the 
early PN phase, if the fast wind speed is moderate, v 2 ~ 400 — 600 km s _1 , and the mass loss 
rate is a fewxlO~ 7 M yr -1 . More generally, the X-ray emission from PNs can be accounted 
for by material ejected at speeds of v 2 ~ 400 — 600 km s^ 1 , whether from the central star 
during the post-AGB or from a companion blowing jets at the very late AGB phase (Soker 
& Kastner 2003). This fast flowing gas hits the slow wind, and passes through a shock 
wave. The post-shock gas is the X-ray emitter. In the next section we study more realistic, 
evolving fast winds in an attempt to further constrain the fast wind properties, which can 
produce observed X-ray emission of PNs. The bipolar X-ray morphology of several observed 
PNs, which indicates an important role of jets rather than a spherical fast wind, cannot be 
explained by the flows studied in this paper, and will be studied with a 2D numerical code 
in a future paper. 

A comment on the mass loss rate and the alternative heat conduction model is needed 
here. Observations of mass loss rates during the PN phase show very low mass loss rates, 
M 2 < 1O~ 8 M yr -1 , and high speeds, v 2 > 1000 km s _1 , (Cerruti-Sola & Perinotto 1989; 
Perinotto et al. 1989). As shown by ASB06, such winds cannot account for the X-ray 
emission as the luminosity will be too low and temperature too high. This led to the idea 
of a heat conduction front between the hot bubble and the cold slow wind gas (Soker 1994). 
The density and temperature in the heat conduction front are intermediate, between those 
of the hot bubble and the cool slow wind gas. Heat conduction can make a significant effect 
and is required to explain the X-ray emission only if after the AGB wind there is a sharp 
transition to the PN-phase wind. Indeed, Zhekov & Perinotto (1996) in calculating the 
effect of heat conduction used a fast wind mass loss rate of M 2 = 7 x 1O~ 9 M yr -1 for a 
fast wind speed of v 2 ~ 500 km s -1 . Namely, the heat conduction model is based on the 
assumption that the relatively high mass loss rate ceases before the post-AGB wind speed 
reaches v 2 ~ 400 km s^ 1 . This type of wind implies a relatively long post-AGB evolution. 

A different scenario would be a dense post-AGB wind, i.e., one in which the mass loss 
rate stays relatively high until later times, at least until the fast wind speed increases to 
v 2 ~ 600 km s" 1 (Soker & Kastner 2003; ASB06). These are the type of winds considered 



- 11 - 



32.6r 



32.4 



i 



w 32.2 

CD 

X 32 

x 



o 31.8 



31.6 



31 



[1.75] 
(400) 



6.2 



[1.17] 

(600) 




6.3 



6.4 



6.5 



6.6 



6.7 



"og(T ) (K) 



Fig. 4. — Comparison of PN evolution tracks in the L x —T x plane between the full simulations 
(+ signs, top to bottom from 40 yrs to 800 yrs, every 40 yrs) and the self similar solutions 
(* signs at 40, 120, 200, 400, and 800 yrs). Radiative cooling was included in the self 
similar calculations by disregarding wind segments with r coo i < t for the X-ray emission. 
Equal-time data points are connected by lines to facilitate the comparison. The three cases 
considered differ by their [M 2 /1O~ 7 M yr _1 ] and (v 2 / km s _1 ) values as indicated. Mi = 
7 x 1O~ 6 M yr _1 and the initial slow wind temperature is 10 K in all runs. Circles (o) in 
the middle run mark the self similar results in which wind segments with r coo i < 3t were 
disregarded. 
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here in Fig. HI Heat conduction, even if it exists here, will not change much, as the density 
of the hot bubble is already high enough for radiative cooling to be important, as seen in 
Fig. [3J This can be demonstrated as follows. At early times the pressure in the hot bubble is 
larger and temperature lower that that assumed in the heat conduction model (Soker 2004). 
Substituting a pressure of nT = 3 x 10 8 cm -3 K and a temperature of Th = 2.5 x 10 6 K, as 
appropriate near the contact discontinuity (see Fig. [3]), in equation (6) of Soker (1994), we 
find that radiative losses become important already at an age of t ~ 500 yr. In the case of 
a very hot bubble (T h = 10 7 K), as assumed by Soker (1994), radiative losses dominate at 
much later times, t > 10, 000 yr. 

There are several reasons to prefer the dense post-AGB wind scenario: (1) Even weak 
magnetic fields, which are likely to exist to some level in all post-AGB stars, will inhibit 
almost completely heat conduction. It is still possible that mixing occurs there (Chu et 
al. 1997), e.g., via turbulence (see the 2D simulations of SS06). (2) A large fraction of the 
X-ray bright PNs are very young (Soker & Kastner 2003), showing that the transition from 
the AGB to the PN phase is short. This implies a high mass loss rate during the post-AGB 
phase. (3) The strong X-ray emitting PN BD +30°3639 (PN G064.7+05.0) currently has a 
wind speed of ~ 700 km s _1 (Leuenhagen et al. 1996) as assumed in the dense fast wind 
model, and much slower than required in the heat conduction model. However, we caution 
that it is possible that the X-ray emission in this PN results from a CFW (jets), and not 
from the central fast wind. To summarize, at this stage both mechanisms in which the X-ray 
emitting gas is blown by the central star should be examined, as we cannot rule out that 
in some PNs the X-ray emission source is a heat conduction front (Soker 1994; Zhekov & 
Perinotto 1996; Steffen et al. 2005). However, at this point we favor the denser post-AGB 
model, which is the subject of this paper. 



4. TIME VARYING FAST WIND 

4.1. General Considerations 

The conclusion of the previous section, which strengthens previous results (Soker & 
Kastner 2003, ASB06, and SS06) is that the fast wind segments that account for the X- 
ray emission should have speeds of vi ~ 400 — 600 km s -1 , and the mass loss rate should 
be M2 ~ 10~ 7 M Q yr -1 . Thus, within a few hundred years the mass of the relevant X-ray 
emitting plasma would be a few times ~ 1O _5 M . The fast wind, however, evolves with 
time. In this section we assume that the X-ray emitting gas in PNs comes from a spherically 
symmetric central stellar wind, and use the X-ray properties of PNs to learn about and 
constrain the evolution of the fast wind. The fast wind evolution depends on the mass of 
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the central star (e.g., Villaver et al. 2002), and possibly on an interaction with a companion. 
There are many uncertainties, and several different functions for the dependence of the fast 
wind mass loss rate and velocity on time have been proposed (e.g., Mellema 1994; Villaver 
et al. 2002; Perinotto et al. 2004; Garcia-Segura et al. 2006). To demonstrate the basic 
properties of an evolving fast wind, we consider several types of fast wind evolution. 

Since slow winds can not produce X-rays, we start our simulations from when the fast 
wind speed is already v 2 — 200 — 300 km s -1 , which is after the PN stage has already started 
(for example, we assume minimum temperature of 10 4 K due to ionization by the central 
star). When comparing to observations, the real dynamical age of the PN, thus, is larger 
than the time (age) of the simulations by several hundred years. Our goal is to conduct a 
general study to find the influence of the different parameters on the X-ray emission, rather 
than to fit each PN with an appropriate mass loss evolution. 



4.2. Different Simulated Cases 

The strategy of varying the wind speed is to take an initial velocity of V20 = 300 km s~ 1 
and to allow it to increase gradually to currently observed speeds. The rise in velocity is 
assumed to take place over a typical time r v = 1000 yr. For simplicity, we start by assuming 
that the rise in velocity is linear. 

V2 = V 20 (l + . (1) 

Our different cases are summarized in Table 1. 

We first consider two versions of wind evolution. 



Table 1: Cases Calculated in section 4.2 



Run 


Mi 


M 2 


Vi 


v 2 






M yr- 1 


M yr" 1 


km s -1 


km s _1 


yr 


A 


10" 5 


3 x lO-^oM) 


10 


eq.l 


1000 


B.l 


10~ 5 


3 x lO-^oM) 2 


10 


eq.l 


1000 


B.2 


10~ 4 


3 x 10- 7 (v 20 /v 2 ) 2 


10 


eq.l 


1000 


B.3 


10" 4 


3 x 10- 7 (v 20 /v 2 ) 2 


10 


eq.l 


500 


B.4 


10- 4 


io-> 2 oM) 2 


10 


eq.l 


1000 


C 


10- 4 


From P04 


10 


From P04 





Notes: (1) t v is defined in equation (1). (2) P04 refers to Perinotto et al. (2004). 
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(A) Constant momentum deposition rate P2 = M 2 v 2 , i.e. decreasing M 2 = 3x 10~ 7 (v 2 o/v 2 )M & yr 
The evolution of the X-ray emission in the L x —T x plane for this case with Mi = 10 _5 M Q yr _1 

(as in all cases here Vi = 10 km s _1 ) is shown in Fig. E]by the '+' symbols. The results are 
shown from t = 250 yr to t = 5000 yr, with marks given every thousand years. It can be seen 
that T x increases with the velocity, while L x decreases since M2 decreases. The total mass 
that is lost to the fast wind over 2000 yr is AM w2 = 3.3 x 10~ 4 M Q . This is approximately 
the mass available for the wind. Therefore, one could expect the wind in this model to cease 
after ~ 2000 yr. 

(B) Constant kinetic power (luminosity) E 2 = (l/2)M 2 v 2 , which together with equa- 
tion (JTJ implies a decreasing mass loss rate: M 2 = M 20 (v 20 /v 2 ) 2 . The total mass lost in the 
fast wind in this case is AM w2 = M 20 t v . We consider several cases with constant kinetic 
power. 

(B.l) First, we take the typical parameter values: M 2 = 3 x 10~ 7 (v 2 o/v2) 2 M Q yr -1 , 
t v = 1000 yr, V20 = 300 km s" 1 and Mi = 10 M© yr -1 . The simulation results are 
represented in Fig. by circles. The temperature T x can be seen to decrease faster in this 
case than in the constant momentum deposition case (A), because here M 2 oc v^ 2 and not 
M 2 oc v 2 \ 

(B.2) Next, we change the mass loss rate in the slow wind to Mi = 10 -4 M Q yr -1 ; all 
other parameters and time scales remain as in run B.l. These results are represented in Fig. 
[S] by asterisks. The denser slow wind confines the hot bubble to a smaller volume, hence 
leading to higher X-ray luminosity. 

(B.3) Next, we want to test a faster evolving wind with a typical time scale of r v = 500 yr 
instead of r v = 1000 yr (equation fl!])). All other parameters are as in (B.2), e.g., the same 
constant wind kinetic power and same slow wind parameters. These results are shown in Fig. 
[5] by star symbols. The faster evolving wind results in a more rapid increase in temperature 
at all times. At early times this leads to higher X-ray luminosity because the post-shock fast 
wind is hotter, emitting stronger in the X-ray band relative to run B.2. At later times the 
luminosity is lower than in run B.2 (for equal times) because the mass loss rate is lower in 
run B.3, and hence the density inside the hot bubble is lower. 

(B.4) The diamonds symbols mark a case similar to B.3, but with denser fast wind: 
M 2 = 10- 6 (v 20 /v 2 ) 2 M & yr- 1 instead of M 2 = 3 x 10- 7 (v 20 /v 2 ) 2 M & yr -1 . As expected, the 
denser bubble has a higher luminosity. 

(C) Finally, we wish to test the wind evolution proposed by Perinotto et al. (2004; 
also Schonberner et al. 2005a,b). Taking v 2 (t) and M 2 (t) from Fig. 1 (Schonberner et al. 
2005a,b) we obtain the results represented in Fig. [5] by the dots. Here v 2 (t) is very low until 
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iog(T x ) (K) 



Fig. 5. — Evolution of X-ray emission from t = 250 yr to t = 5000 yr for the evolving- 
wind cases described in §4.21 Symbols are plotted every 1000 years. Plus signs show the 
constant momentum deposition case (A). Constant luminosity cases B.l - B.4 are represented 
respectively by circles, asterisks, stars, and diamonds. Dots show the evolution for the wind 
according to Perinotto et al. (2004) for which the t = 1000 yr mark is beyond the scale 
of the graph to the left. Triangles represent observed PNs which are listed in Table 2. 
All simulations at late times predict bright, detectable, high-Tx PNs that are not observed 
implying that the wind must evolve (i.e. M 2 {t) must decrease) faster than in these models. 
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t = 250 yr, so we do not see much X-ray until after 1000 yr, hence the 1000 yr mark is not 
on this plot. 

4.3. Comparison with Observations and Discussion 

The various simulations are compared in Fig. [5] to seven observed PNs whose details are 
given in Table 2. It can be seen that a linear wind-speed evolution could generally explain 
the bright PNs during the early phase of t < 1000 yr. For the brightest PNs, the slow 
wind needs to be dense, M\ ~ 10~ 4 M© yr -1 as assumed in runs B.2 - B.4. Also, the high 
temperatures of logTx > 6.5 produced in all of the simulations at late times occur when the 
luminosity is decreasing, but still at a detectable level. However, to date no PN has been 
detected at these high X-ray temperatures. If there are PNs with higher temperatures, they 
must be as faint as the faintest PNs in Table 2 or fainter. This implies that the fast winds 
must evolve faster than assumed in our simulations, so that by the time the temperature 
rises beyond logTx = 6.5, M 2 is too low to produce significant X-rays. 

SS06 found a similar result, although they examined only one case of fast wind evo- 
lution. They conclude that the shocked fast wind cannot account for X-ray properties of 
BD +30°3639 and NGC 40. We disagree with SS06 on that matter, as a more rapidly 
evolving fast wind can in fact account for the X-ray properties. In §4.41 we explore this 
possibility. 

The fast wind evolution as used by Perinotto et al. (2004) results in low luminosities 
inconsistent with six out of the seven observed PNs. The lowest-Lx PN, NGC 40, might 
be explained by this wind. However, the low temperature observed for this target T x ~ 
1.5 x 10 6 K occurs in the Perinotto et al. (2004) wind at ~ 1500 yr, while NGC 40 is 
believed to be 5000 yrs old. 

According to Perinotto et al. (2004), the fast wind speed rises from 200 km s _1 to 
700 km s~ x in 1300 yr. This appears to be too slow. For example, the typical PN BD+30°3639 
has a dynamical age of only ~ 800 yr (Li et al. 2002), but its fast wind speed is already 
~ 700 km s _1 (Leuenhagen et al. 1996), and its mass loss rate is much higher than that 
used by Perinotto et al. (2004) This suggests that the post-AGB evolution of strong X-ray 
emitting PNs is faster than what single star models (Perinotto et al. 2004) predict, but it 
is also possible that the X-ray emission in BD+30°3639 and other similar PNs result from a 
CFW (jets), and not from the central fast wind of a single star. 

There are more reasons to believe rapid pre-PN (post-AGB) evolution. Recent studies 
suggest that the majority of observed PNs are descendant of binary systems (Moe & De 
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Marco 2006; Soker & Subag 2005; Soker 2006). The stellar companion not only is expected 
to shape the AGB wind such that the descendant PN be axi-symmetrical rather than spher- 
ical, but the companion is also likely to disturb the AGB progenitor such that the mass loss 
rate will be much higher (Soker & Subag 2005; Soker 2006). The disturbed AGB envelope 
will affect the post-AGB phase, so that the post-AGB phase will be shorter than that ex- 
pected from single star evolution (Soker 2006). This consideration is another justification for 
invoking a rapidly evolving fast wind ( §4.41) as eventually the mass that is left on the central 
star is very low, and the mass loss rate must decrease fast, reducing the kinetic power in the 
wind. 



4.4. Modified Late-Phase Evolution 

Given the inappropriateness of the previous models at late times, in this section we 
consider a different late-stage evolution scenario in which at some point the mass loss rate 
starts decreasing more rapidly than M 2 ~ v^ 2 ~ t~ 2 . We take the fast wind velocity to 
continue increasing linearly according to equation [H but at time r s the wind's kinetic power 
E 2 is no longer constant, and it starts to decrease linearly: 

E 2 f 1 < t < t s 

E20 [l-{t-T s )r^ r s <t, 

where Td is the typical decay time. We also set a minimum mass loss rate M2 m i n . 

In the first run, we take case B.2, namely M x = 1O~ 4 M yr -1 and t v = 1000 yr at 
early times and for late times using equation [2j we take r s = 500 yr, and = 5000 yr, and 
M 2min = 1O~ 9 M yr" 1 . In the following runs, we vary T4. The results are shown in Fig. [6j 
where we also show one run with M2 m i n = 1O _1O M0 yr -1 and one with M2 m i n = 1O 8 M yr -1 , 
both of them for = 1000 yr. 

The results of Fig. 6 can be interpreted as follows. As the supply of gas decreases, 
the hot bubble density decreases because of expansion, and thus Lx decreases. The energy 
budget of the hot bubble then is dominated by adiabatic cooling. This explains why the 
temperatures are lower than in the previous section, although the evolution of the wind 
velocity v 2 (t) remained unchanged. This type of flow might result from very rapid post- 
AGB evolution, or if the confinement of the hot bubble by the cold shell is not perfect; For 
example, the two jets in NGC 40 (Meaburn et al. 1996) through which hot gas can escape, 
and thus reduce the X-ray temperature and luminosity. 

In Fig. [7] we connect the different PNs to the nearest plot having its age equal to 
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Fig. 6. — Evolution of the X-ray emission for a rapidly decaying fast wind as in equation [2j 
The value of is marked near each run. In all runs Mi = 1O _4 M yr -1 , t v = 1000 yr, 
and t s = 500 yr. M 2m i n = 1O~ 9 M yr -1 in all cases except the one marked 10~ 10 and the 
one marked 10~ 8 . Symbols mark 1000 yr intervals, up to 4000 yr in the run marked with 
diamonds, and up to 5000 yr in all other runs. All plots start at 250 yr. 
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the dynamical age of the PN. The dynamical age of the PNs are marked on the figure. We 
removed three runs and added one run where we started the fast wind with v 2 o = 200 km s -1 , 
and took r s = 1000 yr, r d = 500 yr, M 20 = 10" 6 M Q yr _1 , and M x = 1O" 4 M yr" 1 (star 
symbols). This figure shows that in addition to the luminosity and the temperature of the 
X-ray emitting gas, we can account for the age of PNs as well. 

Since we start the simulations at the PN stage, the PN time t in the simulation is 
somewhat smaller than the actual age of the PN, by ~ few x 100 yr. Therefore, in Fig. [71 
we connected the observed PNs to simulated t values that are somewhat smaller than the 
observed age of the PN. For example, the upper left PN whose estimated age is 1000 yr 
is connected by a line to a t — 750 yr simulation data point. We mention again that our 
goal is to find the general constraints on the properties of the fast wind, if it is the source 
of the X-ray emission, rather than to fit individual PNs. Changing a little the parameters 
used here, we can fit each PN individually. In addition, the composition of the fast wind 
also affects the X-ray properties (SS06); we use constant (solar) composition for the cooling 
function. 

We finish this section by showing the evolution with time of the most important physical 
quantities. In Fig. [8] we show the evolution with time of the quantities in the run with M\ = 
1O" 4 M yr" 1 , t v = 1000 yr, t s = 500 yr, r d = 1000 yr, and M 2min = 1O- 1O M yr -1 . The 
total mass lost in the fast wind (from t = in our calculation, when v 2 = 300 km s^ 1 ) at 500, 
1000, and 1500 years is AM 2 = 1O~ 4 M , 1.387x 1O~ 4 M , and 1.4675 x 1O~ 4 M , respectively. 
Initially, the reverse (inner) shock moves outward, as does the contact discontinuity. At 
t = 1250 yr it starts to move inward because the ram pressure of the fast wind decreases by 
a large factor. At t = 1750 yr it starts to move outward again. 

5. SUMMARY 

Our main finding is that if the source of the extended X-ray emitting gas in PNs is the 
fast wind blown by the central star, then the fast wind must be rapidly evolving. Otherwise 
the luminosity would be too low at early times, and the fast wind would produce hot gas at 
later times that is not observed. Our simulations require that the mass loss rate of the fast 
wind be relatively high, M 2 = 1 — 3 x lO~ 7 A<f yr -1 , at early times when the wind speed 
is ~ 300 — 700 km s _1 , and then it needs to rapidly decline to (Fig. [H]). This weakening 
of the wind occurs approximately 500 — 1500 yr after the fast wind has started to blow, or 
~ 1000 — 3000 yr after the star has left the AGB. Our results are contradictory to those 
of SS06, who examined only one case of fast wind evolution, and could not find a match 
between their model and the X-ray properties of BD +30°3639 and NGC 40. 
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Fig. 7. — As Fig. [6j with three runs omitted. The dynamical age of PNs (in years)are 
marked, and each PN is connected to a wind age equals to its dynamical age. The stars 
marked a new run with parameters as written on the upper left corner. 
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Fig. 8. — From top to bottom, the evolution of the fast wind speed in km s" 1 , the mass 
loss rate in 1O _7 M yr^ 1 , the logarithm of the ratio of kinetic power to its initial value, the 
X-ray luminosity in 10 -32 erg s _1 and temperature in 10 6 K, for the run marked by squares 
in Figs. Eland [7]). The first point of each run shown on Figs. Eland [7] is at an age of 250 
years, as marked here by the dashed vertical line. Note that because of radiative cooling of 
cooler gas, T x increases at late times. 
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A rapid fast- wind evolution is compatible with a recent claim for rapid post-AGB evolu- 
tion caused by binary interaction (Soker 2006). The high mass loss rate at the early PN phase 
is a continuation of a high mass loss rate during the final AGB phase. Such a high mass loss 
rate may lead to dust formation and strong IR emission. Montez (2006; R. Montez, private 
communication) suggested a correlation between IR emission and X-ray emission. Not all 
PNs will have such rapid fast-wind evolution, hence not all PNs will have detectable X-ray 
emission. Other PNs, those with slowly evolving fast winds are not expected to produce 
detectable X-ray emission as adiabatic cooling overcomes shock heating. Note that a slowly 
evolving fast wind requires a low mass loss rate. The reason is that during the post-AGB 
phase the mass in the envelope is very low, and the star cannot lose mass at a high rate for 
a long time. The low mass lose rate implies a weak X-ray emission, as the case C in Figure 

m 

We should stress here that we are not claiming that all PNs owe their X-ray emission to 
the central star fast wind. Definitely, in some cases jets are the source of the X-ray emitting 
gas. If observations indicate that the fast wind evolution in PNs is much slower than required 
by our models, (e.g., Perinotto 2004), we will have to agree with SS06 and conclude that 
a central star wind can not explain the observed X-rays. Jets blown by a companion may 
be an alternative explanation for the X-ray emission. In this case we also expect that only 
a fraction of all PNs will have detectable X-ray emission. Heat conduction can also play a 
role, although our models show that heat conduction is not required if the fast wind evolves 
fast enough. 

As both rapid fast-wind evolution and jets require the presence of a close companion 
(in the case of jets a stellar companion, in the case of rapid evolution a brown dwarf or even 
a massive planet might do the job as well), we propose that the extended X-ray emission in 
PNs is tightly connected to the presence of a close binary companion. 

We thank Raghvendra Sahai, Matthias Stute, Martin Guerrero, Joel Kastner, and an 
anonymous referee for useful comments. This research was supported in part by the Israel 
Science Foundation, grant 28/03, and by the Asher Fund for Space Research at the Technion. 
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Table 2: X-ray properties of Planetary Nebulae 



# 


PN 


L x 


T x 


Dynamical Age 






10 32 erg s _1 


10 6 K 


yr 


1 


NGC 7027(PN G084.9-03.4) 


1.3 


3 


600 


2 


BD +30 3639(PN G064.7+05.0) 


1.6 


3 


700 


3 


NGC 7026(PN G096.4+29.9) 


4.5 


1.1 


< 1000 


4 


NGC 6543(PN G096.4+29.9) 


1.0 


1.7 


1000 


5 


NGC 7009(PN G037. 7-34.5) 


0.3 


1.8 


1700 


6 


NGC 2392 (PN G197.8+17.3) 


0.26 


2 


1800 


7 


NGC 40 (PN G120.0+09.8) 


0.02 


1.5 


5000 



The parameters of the first three PNs and NGC 7009 are summarized by Soker & Kastner (2003). The data 
for NGC 2392 are from Guerrero at al. (2005), for NGC 40 from Montez et al. (2005), and for NGC 7026 
from Grucndl ct al. (2006). 



